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Late internal shock model for bright X-ray flares in 
Gamma-ray Burst afterglows and GRB 011121 
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ABSTRACT 

We explore two possible models which might give rise to bright X-ray flares in GRBs 
afterglows. One is an external forward-reverse shock model, in which the shock param- 
eters of forward/reverse shocks are taken to be quite different. The other is a so called 
"late internal shock model" , which requires a refreshed unsteady relativistic outflow 
generated after the prompt 7— ray emission. In the forward-reverse shock model, after 
the time t x at which the RS crosses the ejecta, the flux declines more slowly than 
(ie/i x )~^ 2+/5 \ where i® denotes the observer's time and (3 is the spectral index of 
the X-ray emission. In the "late internal shock model" , decaying slopes much steeper 
than (i©/i c ,©)~^ 2+ ^' ) are possible if the central engine shuts down after t 6iffi and the 
observed variability timescale of the X-ray flare is much shorter than i e ,®. 

The sharp decline of the X-ray flares detected in GRB 011121, XRF 050406, GRB 
050502b, and GRB 050730 rules out the external forward-reverse shock model directly 
and favors the "late internal shock model". These X-ray flares could thus hint that 
the central engine operates again and a new unsteady relativistic outflow is generated 
just a few minutes after the intrinsic hard burst. 
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1 INTRODUCTION 

GRB 011121 was simultaneously detected by BeppoSAX 
GRBM and WFC (Piro 2001), and the fluence in the 2- 
700 keV range corresponds to an isotropic energy of 2.8 x 
10 52 ergs at the redshift of z = 0.36 (Infante et al. 2001). This 
burst was born in a stellar wind (Price et al. 2002; Greiner et 
al. 2003) and a supernova bump was detected in the late op- 
tical afterglow (Bloom et al. 2002; Garnavich et al. 2003). Its 
very early X-ray light curve, which has not been published 
until quite recently, is characterized by the presence of two 
flares (Prio 2005, hereafter P05). In the first one, which is 
also the strongest of the two, the observed flux F rises and 
decays very steeply: F oc t^ for 239 s < t® < 270 s and 
F oc t^ 7 for 270 s < % < 400 s, where t® is the observer's 
time 1 . Such a peculiar flare in the early X-ray lightcurve of 



* E-mail: yzfan@pmo.ac.cn(YZF); dmwei@pmo.ac.cn(DMW) 
1 During our revision, bright X-ray flares peaking a few minutes 
after XRF 050406, GRB 050502b and GRB 050730 have been 
reported (Burrows et al. 2005; Starling et al. 2005). Sharp rise 
and fall are also evident in these events. 



GRBs has not been predicted before. P05 suggested that the 
X-ray flare represents the beginning of the afterglow. 

In this Letter we explore two alternative models which 
might give rise to very early X-ray flare in GRB afterglows 
(§2): a forward-reverse shock model (§2.1) and a "late in- 
ternal shock model" (§2.2). We compare the available data 
with the predictions of those models in §3 and summarize 
our results in §4, with some discussion. 



2 POSSIBLE MODELS 

2.1 The external forward-reverse shock model 

The external forward-reverse shock model has been widely 
accepted on interpreting the early IR/optical flashes of GRB 
990123, GRB 021211 and GRB 041219a (For observations, 
see: Akerlof et al. 1999; Fox et al. 2003; Li et al. 2003; 
Blake et al. 2005. For theoretical modeling, see: Sari & Piran 
1999; Meszaros & Rees 1999; Wei 2003; Kumar & Panaitescu 
2003; Fan, Zhang & Wei 2005). Synchrotron radiation from 
the reverse shock (RS) and the forward shock (FS) usually 
peaks in the infrared-to-optical and ultraviolet-to-soft X- 
ray bands, respectively. Thus, the RS emission component 
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is not dominant in the X-ray band. The synchrotron self- 
Compoton (SSC) scattering effect of the RS radiation has 
also been considered by different authors, but no strong X- 
ray emission is found to be expected (Wang, Dai & Lu 2001) 
except in some carefully balanced conditions (Kobayashi et 
al. 2005). 

In most of previous works, the fractions of FS energy 
given to electrons, e B , and to magnetic field, e B , were as- 
sumed to be the same as the corresponding fractions in the 
RS. However, this may not necessary be the case. Fan et al. 
(2002) performed a detailed fit to the optical flash of GRB 
990123 data and obtained e r c = 4.7el and e B = 400e B , where 
the superscripts "r" and "f" represent RS and FS, respec- 
tively. Similar results were obtained by Zhang, Kobayashi & 
Meszaros (2003), Kumar & Panaitescu (2003), Panaitescu 
& Kumar (2004), McMahon, Kumar & Panaitescu (2004), 
and Fan et al. (2005). In this section, we study the RS/FS 
emission in X-ray band by adopting different shock param- 
eters. We focus on the thin shell case (i.e., the RS is sub- 
relativistic, see Kobayashi [2000]), in which the RS emission 
is well separated from the prompt 7— ray emission. 

ISM model. In the thin shell case, the observer's time 
at which RS crosses the ejecta can be estimated by (e.g., 
Fan et al. 2005) 

tx « 128s( i^ )£ V3 B3n -l/3, ? -8/3 ; (1) 

where Ei SO is the isotropic energy of the outflow, n is the 
typical number density of ISM, and r\ is the initial Lorentz 
factor of the outflow. The convention Q y = Q/10 y has been 
adopted in cgs units throughout the text except for some 
special notations. 

In the standard afterglow model of a fireball interacting 
with a constant density medium (e.g., Sari, Piran & Narayan 
1998), the cooling frequency fc,©, the typical synchrotron 
frequency and the maximum spectral flux -F^ max read: 

< e = 4.4 x 10 17 Hz i^&e^fo V-3 (pX 1 + 

= 4.4 x 10 15 Hz 4o,b 3 4,-3 ee 2 -i*d,-3^(T^), and 

^ f ,ma X = 2.6 mjy £ iS o,534,-3 n o - D L,28.34( i T £ ) ) where 
Cp = 13(p — 2)/[3(p — 1)], p ~ 2.3 is the typical power- 
law distributio n index of th e electrons accelerated by FS, 
Y t ~ [—1 + yl + 4£ f di/e B ]/2 is the Compton parameter, 
x i w min{l, (^/^') (p_2)/2 } (Sari & Esin 2001), and D L is 
the luminosity distance for (fi M , Ha, h) — (0.3, 0.7, 0.71). 
Hereafter t = t©/(l + z), and td is in unit of days. 

Following Zhang et al. (2003), we take ej, = K c e c and 
e B = TZ%e B . Att x , the RS emission satisfies [See also Fan et 
al. (2005), note that a novel effect taken into account here 
is the inverse Compton cooling of the electrons] 

*4, e (*x) = ^b[^ c ( 73 4,x - i)] 2 *4,©(t x )/(r x - 1) 2 (2) 

<©(**) « K B 3 [(l + Y f )/(l + Y r )] 2 ^(t x ), (3) 

f; m »(tx) « ^ B < max (i x ). (4) 

where 734.x ~ (»)/r x + T x /r))/2 is the Lorentz factor of 
the shocked ejecta relative to the initial one, F x is the 
bulk Lorentz factor of the shocked ejecta at tx, Y T ~ 
[— 1 + + 4x r K c el/(K B e B )]/2 is the Compton parameter 
and x 1 w min{l, (i4/i/ c r ) (p - 2)/2 }. 

If both v r c (S and z>£© are below the observed frequency 
f© , the detected flux of RS and FS emission are 



K 9 (t x ) ~ F;, max (t x )K 1 , ffi (t x )]( p - 1 )/ 2 [< ffi (tx)] 1/2 ^ p/2 ,(5) 

Ft 9 (tx) ~ < max (t x )[^,©(tx)] (p - 1)/2 [<©(tx)] 1/2 ^ p/2 .(6) 

Taking p — 2.3 (such a choice is to match the ob- 
served slope of the X-ray flare f) = p/2 — 1.15, see P05), 
T x « r?/2 ~ 100, e Bj _ 3 = 1, = 1, Kb = 10 and 

ft e = 5, we have x i w 1, Y f w (e c /e B ) 1/2 = 10, a; r « 0.6, 
Y" r w 1.2, and (1 + F f )/(1 + y r ) > 1, i.e., we have a larger 
contrast F^, (tx) / (tx) when the inverse Compton effect 
has been taken into account. With equation (7), we have 
Fl (t x ) I 'Fl (t x ) ~ 5, i.e., in the X-ray band, the RS emis- 
sion component is dominant. For £© > tx, the RS emission 
declines as (t®/t x ) ~( 2+p / 2 ^ because of the curvature effect 
(e.g., Kumar & Panaitescu 2000, hereafter KP00) 2 . The FS 
emission declines as t^~ 3p ^ 4 (e.g., Sari et al. 1998), so the 
X-ray flare lasts ~ [FZ 9 (t x )/Fl B (tx )] 4/(1 °- p) t x ~ 300 s. 
Taking z = 0.36, Q y = 1, and v x = 2.4 x 10 17 Hz, we 
have VxFl^ ~ 10~ 9 ergs cm" 2 s" 1 [(1 + z)/1.36]Dl 2 27 . 7 . 
The peak flux of the X-ray flare is ~ u x [F^ + F^ x ] ~ 
5 x 10 -9 ergs cm~ 2 s _1 [(1 + z)/1.36]Dl i 2 . 7 .7, which is con- 
sistent with the observation of GRB 011121 (~ 2.4 x 
10 -9 ergs cm -2 s _1 , see P05). However, the accompany- 
ing optical flash is very bright. With the typical parameters 
taken here, the V band flux is ~ 5 Jy. 

Wind model. GRB 011121 was born in a stellar wind. 
The best fit parameters are p = 2.5, -Eiso, 52 = 2.8, A, ~ 
0.003, el ~ 0.01, and e B ~ 0.5 (P05). It is straightforward 
to show that with proper choice of K e and Kb, at tx, the 
RS emission may be dominant in the soft X-ray band. 

Numerical results. Following Fan et al. (2005), the 
FS-RS emission (in the X-ray band) has been calculated 
numerically. In the ISM case (see Fig. 1(a)), the parameters 
are taken as £ iso ,53 = 1, p = 2.4, 4 = 0.1, e B = 0.001, 
n — 1 cm~ 3 , and the initial width of the outflow A = 6x 10 11 
cm. In the wind case (see Fig. 1(b)), we take the best fit 
parameters presented in P05. 

As shown in Fig. 1(a), in the ISM case, there comes an 
X-ray flare dominated by the RS emission only when both 
Kb and K c are much larger than unity. In the wind case, 
with proper Kb and K c , the RS emission may be dominant 
in the soft X-ray band, too. But there is no flare expected 
since both the FS and RS emission components decrease 
continually even at very early time (see also Zou, Wu & Dai 
2005). So the FS-RS scenario is unable to account for the X- 
ray flare detected in GRB 011121. Moreover, in the general 
framework of accounting for X-ray flares in GRB afterglows, 
the FS-RS model is further disfavored for the fact that, even 
in an ISM scenario with parameters suitable for a RS flare 
arising above the FS emission, the predicted temporal decay 

2 To derive the curvature effect, two assumptions arc made 
(KP00). One is that the Lorentz factor of the outflow is nearly 
a constant. The other is that the observer frequency should be 
above the cooling frequency of the emission. As far as the reverse 
shock emission mentioned here, these assumptions are satisfied. 
So we take (t©/tx )~( 2+p / 2 ) to describe the decline, which has 
been verified by the detailed numerical calculation (Fan, Wei & 
Wang 2004; sec also Fig. 1 of this work). 
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Figure 1. The very early X-ray (u^ = 2.42 X 10 17 Hz) sample 
lightcurves. For the reverse shock emission component, TZ e and 
have been marked in the figure, (a) The ISM case, the param- 
eters are taken as n = 1 cm -3 , -Ei SOj 53 = 1, z = 1, A = 6x 10 11 
cm, 4 = 0.1, 4 = 0.001, r? = 200,' and p = 2.4. (b) The wind 
case, following P05, we take A, = 0.003, 4 = 0.01, t l B = 0.5, 
Aso,53 = 0.28, p = 2.5, and z = 0.36; In addition, we assume 
r] = 200 and A = 3.0 X 10 12 cm. 

is too shallow compared to the steep decay of X-ray flares 
observed so far (see Fig. 1). 

2.2 Late internal shock model 

In the standard fireball model of GRBs, the 7— ray emission 
is powered by internal shocks, whose duration depends on 
the active time of the central engine. However, the variability 
of some GRB afterglows implies that the activity of the GRB 
central engine may last much longer than the duration of 
the prompt emission recorded by 7— ray monitors (e.g., Dai 
& Lu 1998; Granot, Nakar & Piran 2003; Ioka, Kobayashi 
& Zhang 2005). In addition, it has been proposed that the 
Fe line observed in some GRB X-ray afterglows could be 
attributed to a prolonged activity of the central engine (Rees 
& Meszaros 2000; Gao & Wei 2005). 

A possible mechanism for the re-activity of the central 
engine could be as follows. During the accretion phase which 
powers the prompt 7— ray emission, a fraction of the mate- 
rial constituting the massive progenitor could possibly be 
pulled out; the central engine could thus be restarted at late 
times by the fall back of part of this material onto the central 
collapsar remnant (King et al. 2005) . 

Here we assume that the central engine restarts a few 
minutes after the prompt 7— ray emission, powering a new 



unsteady relativistic outflow. We suppose that the Lorentz 
factor of the ejected material can be highly variable, set- 
ting r s ~ 10 and Tf ~ 100, as the typical Lorentz factors 
of the slow and fast shells, respectively. The masses of the 
slow and fast shells are taken as m,{ ~ m a . When an inner 
fast shell catches up with an outer slow shell at a radius 
~ 2Tlc8t<$) I '(1 + 2) (where 8t@ is the observed typical vari- 
ability timescale of the X-ray flares) , internal shocks are gen- 
erated. The Lorentz factor of the merged shell is T m \/r f T s 
(e.g., Piran 1999), and the Lorentz factor of the internal 
shocks can be estimated by T sh « (y/Ti/T e + y/T a /T t )/2. 
We call this re-generated internal shocks as the "late internal 
shocks" . 



2.2.1 Physical parameters 

The internal shock model has been discussed by many 
authors (e.g., Paczyhski & Xu 1994; Rees & Meszaros 
1994; Daigne & Mochkovitch 1998, 2000; Piran 1999; Dai 
& Lu 2002). Generally speaking, to calculate the syn- 
chrotron emission of internal shocks, the following param- 
eters are involved: the outflow luminosity L m , <5t®, Y, r sh , 
and of course the shock parameters e c and eb- For typi- 
cal GRBs, (L m , 5t®, T, r s h, e e , €b) are taken to be ~ 
(10 52 ergs s -1 , 0.001-0.01 s, 100-1000, a few, 0.5, 0.01- 
0.1), respectively (e.g., Dai & Lu 2002). 

The time averaged isotropic luminosity (2-700 keV) 
of the X-ray re-bursting of GRB 011121 is L x ~ 6 x 
10 48 ergs s" 1 (P05). So we normalize our expression by tak- 
ing _L m ~ 5x 10 49 L x ,49(e/0.2) 1 ergs s 1 , where e is the effi- 
ciency factor of the X-ray flare. We take e ~ 0.2, as found in 
typical GRBs/XRFs (Lloyd- Ronning & Zhang 2004). Such 
small Lm (comparing with the GRB case) implies that the 
fallback accretion rate is just ~ 0.001 — 0.01 times that of 
the prompt accretion, if the efficiency factor of converting 
the accretion energy into the kinetic energy of the outflow 
is nearly a constant (MacFadyen, Woosley & Heger 2001). 

The <St® measured in X-ray flares is significantly longer 
than that measured in the intrinsic hard burst (Burrows et 
al. 2005). In this Letter, we take 5t® ~ 10 s. The spectra of 
the X-ray flares detected so far are all nonthermal. The op- 
tical thin condition implies a lower limit on the bulk Lorentz 
factor of the merged shell (e.g., Rees & Meszaros 1994) 

r > 11 ^,4 9 . 7 [1-36/(1 + z)]-^St~f>, (8) 

The typical radius of the late internal shock is Tiint ~ 
2T 2 c8t S) /(l + z) > 5 x 10 13 cm L^ /5 49 7 [1.36/(1 + z)] 3 / 5 ^. 

As we will show in the following section, with L m ~ 
5 x 10 49 ergs s" 1 , 5t® ~ 10 s, T ~ 30, r sh ~ 2, e c ~ 0.5, and 
£b ~ 0.1, most of the internal shocks energy is emitted in 
the soft X-ray band and the predicted flux matches the ob- 
servation of the X-ray flare in GRB 011121. Alternatively, as 
shown in Barraud et al. (2005), with e ~ 0.01 (correspond- 
ingly, L m ~ 10 51 - 10 52 ergs s" 1 and r sh ~ 1.06), St 9 ~ 
a few seconds, and T ~ a few hundreds, X-ray dominated 
luminosity is expected, too (Please note that their calcula- 
tion of the internal shocks emission is somewhat different 
from ours; see Barraud et al. 2005 for detail). Therefore, we 
believe that with proper parameters, X-ray flares do appear 
in the late internal shocks scenario. 
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2.2.2 The synchrotron radiation of the "late internal 
shocks" 

Following Dai & Lu (2002), the comoving number den- 
sity of the unshocked outflow is estimated by n c w 
Lm/(4ivT 2 Rf nt mpC 3 ), where m p is the rest mass of proton. 
The thermal energy density of the shocked material is cal- 
culated by e w 4r s h(r s h — l)n e m p c 2 (Blandford & McKee 
1977). The intensity of the generated magnetic field is es- 
timated by B « (87re B e) 1/2 w 6 x 10 3 G [r sh (r sh - 

i)/2] 1/2 ii / , 2 49 . 7 r r.^r nt 1 , 14 . 5 - 

As usual, we assume that in the shock front, the accel- 
erated electrons distribute as dn e /d7 e oc 7<T P for 7 C > 7 e ,m, 
where 7 c>m = e c (r sh — l)[(p — 2)m p ]/[(p — l)m c ] is the min- 
imum Lorentz factor of the shocked electrons (Sari et al. 
1998), and m c is the rest mass of electron. In this section, 
we take p=2.5. The observed typical frequency of the syn- 
chrotron radiation reads 



= 7e m9erB/[2(l + 2)7rm e c] 



2.7 x 10 1D Hz <_o. 3 4(-i(rsh - l) s/2 (r sh /2) 



,5/2 ( 



,1/2 



-1/2 



7ri.5^©,i, 



(9) 



i.e., most of the shock energy is emitted in the soft X-ray 
band, where g e is the charge of electron. 

The cooling Lorentz factor is estimated by (e.g., Sari et 
al. 1998) 7c, c w 7.7 x 10 8 (1 + z)/(TB 2 8t (S ), and the corre- 
sponding cooling frequency reads 



^c,® = 7 0jC g e rS/[2(l + z)7vm c c] 

~ 10 10 Hz [1.36/(1 + 2 )] 2 e B 3 _ /2 r?. 5 

[r sh (r sh -i)/2]- 3/2 L m % 2 7 5t ffi ,i. 



(10) 



The synchrotron self-absorption frequency is estimated by 
(Li & Song 2004) 



z/ a ,© ~ 10 iJ Hz [1.36/(l + z)] 3// L 



3/7 T 2/7 



r- 5/7 sr 4/7 B 1/7 (id 

m.49.7 1 1.5 0L ®,1 D 3 -K 11 ) 

The maximum spectral flux of the synchrotron radiation 
is (e.g., Wijers & Galama 1999) F max « 3V3$ P (1 + 
z)N c m c c 2 a T TB/(32iT 2 q e D'l), where iV e = L m St/[(l + 
z)rm p c 2 ] = 8 x 10 51 [1.36/(1 + z)]L mA9 . 7 r~lSt 9:1 is the 
number of electrons involved in the emission. $ p is a func- 
tion of p, for p = 2.5, $ p w 0.6 (Wijers & Galama 1999). 
For u c>(S < i/ a ,e < ^ x < I'm,©, the predicted flux is (e.g., Sari 
et al. 1998) 



F — F 



(fm,e/^c,e) (^x/fm,©) 



-p/2 



« 2.5 mJy [i/ x /(2.42 x 10 17 Hz)]- p/2 e p - o^I 2 ? 74 

(r sh /2)< p - 2 )/ 4 (r sh - i)< 5p - fl >/ 4 L£+T 
r 2 7<7 )/2 £>L 2 2 7.7- (12) 

Taking Q y = 1 and u x = 2.42 x 10 17 Hz, with equation (12) 
we have F Vyl w 2.5 mJy, which matches the observation of 
GRB 011121 (~ 1 mJy). The V band flux can be estimated 
as (z^c,© < v w < z/ a ,©) F„ v ~ ^max^l^©^ 72 ~ 40 mJy. 

What happens after the "late internal shocks"? Surely, 
the refreshed relativistic outflow will catch up with the ini- 
tial outflow when the latter has swept a large amount of 
material and got decelerated. That energy injection would 
give rise to a flattening (e.g., Rees & Meszaros 1998) or 
re-brightening signature (e.g., Panaitescu, Meszaros & Rees 




Figure 2. The X-ray lightcurve of one flare consisting of 
ten pulses (the solid line, for illustration), each takes the 



profile F v 



*eje)/<5%] for *eje < *a 



and F„ x oc [(t© - t cjc )/<5t©]" 



+ 

5t©. In these pulses (j = 1 — 10), the peak of F VyL ^ arc 
taken to be (0.3, 0.8, 0.6, 0.9, 0.1, 0.7, 0.5, 0.3, 0.7, 1.0), 
respectively (in arbitrary units); 5t©,j are taken to be 
(10, 15, 5, 10, 13, 5, 9, 11, 5, 60) s, respectively; t c jc,i are taken 
to be (100, 110, 125, 130, 140, 153, 158, 169, 185, 190) s, respec- 



tively. t e 



*eje,lO + <5%,10 = 250 s. The dashed line represents 



(W*e,ffi)" 



1998; Kumar & Piran 2000; Zhang & Meszaros 2002), 
which could potentially account for the late re-brightening 
of XRF 050406, the late X-ray re-bursting detected in GRB 
050502b and the second/weaker X-ray bump observed in 
GRB 011121. However, the detailed lightcurve modeling is 
beyond the scope of this Letter. 



2.2.3 The decline behavior of the flare 

For observer's frequencies above the cooling frequency, 
the curvature effect dominates the temporal behavior of 
the observed flux after the central engine shuts down 
at t , ffi (PK00). The flux declines as XX*,4(*(B - 



ieje.i) / <5t©, 



l-(2+ Pi /2) 



(PK00), where i represents the ith 



pulse, t e j e ,i and t©,i are the ejection time and the variabil- 
ity timescale of the iih pulse, respectively. Such a decline is 
much steeper than (t©/£o,©)^ 2+p/ ' 2 ' ) for i e ,® ^> max{5t© ji }. 
For example, as shown in Fig. 2, the decline of the flare 
is dominated by the curvature effect of the last long pulse 
with <5f© ~ 0.24i Cj ©. A crude power-law fit to the decline 
yields F Vx oc (i©/ic,©)~ s , which is steep enough to match 
the sharpest decline detected so far. In reality, the central 
engine does not turn off abruptly. The dimmer and dimmer 
emission powered by the weaker and weaker "late internal 
shocks" may dominate over the curvature effect of the early 
pulses, resulting in a shallower decay. 
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3 DECLINE BEHAVIOR OF THE X-RAY 
FLARE: CONSTRAINT ON THE MODEL 

Early X-ray flares have been well detected in GRB 011121, 
XRF 050406, GRB 050502b and GRB 050730 (P05; Burrows 
et al. 2005; Staring et al. 2005). The rise and fall of the first 
flare (also the dominant one) in GRB 011121 are both very 
steep. Similar temporal behavior is evident in other events. 
The sharp decline of these flares imposes a robust constraint 
on the model, as shown below. 

The X-ray flare detected in GRB 011121 appears at 
th,(B = 239 s and peaks at £ Pi = 270 s. The burst is believed 
to be born in a weak stellar wind. As shown in Fig. 1(b), 
no flare is expected in the FS-RS model. The FS-RS shock 
model is further disfavored by its shallow decline. In the late 
internal shock model, the decline of the flare can be steep 
enough to account for the observation (see Fig. 2 for illustra- 
tion). Moreover, as shown in §2.2.2, with proper parameters 
the observed flux can be well reproduced. So the "late inter- 
nal shock model" is favored. We would like to point out that 
the fall of the X-ray flare detected in GRB 011121 is still at- 
tributed to the late internal shocks rather than the curvature 
effect. The reason is as follows. Since 5t® < (i P ,e — *b,©) = 31 
s, the resulted decline F„ x oc [l + (i®-270)/<5£©]~ 3 15 is much 
steeper than the observation F Vx oc [(t® — 239)/31]~ 1,4 (see 
Fig. 7 of P05) if after t p ,® there are no internal shocks any 
more. 

The X-ray flare detected in XRF 050406 peaks at i Piffl w 
210 s and declines as F oc t^ 5 ' 7 . The X-ray flare detected 
in GRB 050502b peaks at i Pi ® « 650 s and declines as F oc 
t^ r . In the X-ray afterglow lightcurve of GRB 050730, there 
are three X-ray flares (ranging from 200 s to 800 s after 
the trigger of the GRB). A crude fit to the decline of these 
three flares results in F oc t^ 5 or steeper. Obviously, the 
FS-RS scenario is ruled out by the steep observed decays 
and the "late internal shock model" is favored. For the X- 
ray flare detected in GRB 050502b, the late internal shock 
model interpretation is further supported by the sharp spike 
detected in 1.0-10.0 keV band (Burrows et al. 2005). 



4 SUMMARY & DISCUSSION 

In this work, we have explored two possible models which 
might give rise to X-ray flares in GRB afterglows. One is 
the external forward-reverse shock model (the ISM case), in 
which the shock parameters of forward/reverse shocks are 
taken to be quite different. The other is the "late internal 
shock model", which requires that a refreshed unsteady rela- 
tivistic outflow is generated after the prompt 7— ray emission 
(see Ramirez-ruiz, Merloni & Rees [2001] for alternative sce- 
narios) , perhaps due to the fallback accretion onto the central 
collapsar remnant. The refreshed outflow may be character- 
ized by a low outflow luminosity (~ 10 49 ergs s _1 ), a small 
bulk Lorentz factor (~ 30), and a long variability timescale 
(~ 10 s) . In the external forward-reverse shock model, after 
the peak of the reverse shock emission (i Pl © = t x ), the flux 
can not decline more sharply than (teAp.e) _ ^ 2+p ^ 2 - ) (see 
Fig. 1 for illustration). In the "late internal shock model", 
the decline can be much steeper than (t® /i c ,e) _ ' 2+p ^ 2 - ) if the 
central engine shuts down after t Bj $ and the longest variabil- 
ity timescale of the X-ray flare is much shorter than t Cj e (see 
Fig. 2 for illustration). 



For the X-ray flares detected in GRB 011121, XRF 
050406, GRB 050502b and GRB 050730, the external 
forward-reverse shock model is ruled out directly by its shal- 
low temporal decay. For the same reason, other possible ex- 
ternal models (i.e., the model related to the external for- 
ward shock), including the density jump model, the two- 
components jet model, the patch jet model as well as the 
energy injection model are ruled out too (Zhang et al. 2005). 
Thus, the "late internal shock model" is found to be favored. 
In this model, the optical emission may be suppressed due 
to strong synchrotron-self-absorption. But in the ultraviolet 
band, the radiation could be quite strong. Large amount of 
neutral gas would be ionized, as detected in GRB 050502b 
and GRB 050730 (Burrows et al. 2005; Starling et al. 2005). 

Very early X-ray flares are well detected both in long 
GRBs and in XRFs, which strengthens the correlation of 
these two phenomena, though the nature of XRFs is still 
unclear (Barraud et al. 2005 and the references therein). 

Finally, we suggest that the early X-ray light curve of 
some GRBs may be a superposition of the emission powered 
by the long activity of the central engine and the emission 
of the external forward shock. As a consequence, the X-ray 
temporal behavior may be quite different from that of the 
long wavelength emission (UV/Optical ones). This predic- 
tion can be tested by the UVOT and XRT on board Swift 
observatory directly in the near future. 
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face" R(0) = c(t® - i cjc )/[(l + z)(l - VcosO)] (Rees 1966, 
Nature, 211, 468), the magnetic field, the typical emission 
frequency of the electrons, and the number of electrons in- 
volved in the radiation are all functions of 7?, that's why we 
take into account the term (i?/i? cro ) k in calculating S v i. 
Equation (Al) yields 



(1 - Vcos(9i) k (l - Vcos6y {3+p+k) smOdO 



2 + P + k 



(1 - Vcos 9 t ) 



-(2+/3) 



(A2) 



On the other hand, F v{to) oc ^W+ki 1 ~ Vy (2+fJ) . For to < 
A - Vcos9, 



< tj, we have 



(to) I 



^i-(2+/3) 



v 



cjoN-(2+/3) 



(A3) 



which coincides with that presented in KP00 but derived in 
a different way. 3 



APPENDIX A: DERIVATION OF THE 
CURVATURE EFFECT 

Assuming a shell is ejected at i j and moves with a constant 
Lorentz factor T (the corresponding velocity is V, in unit 
of the speed of light), the electrons are shock-accelerated 
at R < Rao and cools rapidly, where R cro is the radius 
after which there is no newly relativistic electrons injected. 
For R > i?cro, the radiation nearly cuts off as long as the 
observer frequency z/q is above the cooling frequency of the 
electrons. For i® > to = t c j c + (1 + z)R CIO /(2T 2 c), the flux 
received from the shell is given by 



f 



, S„>sm9d9 . . 

F ^ (tffi)0C / r3(i-vcosg)3 ' (A1) 



where 9 t satisfies R CIO ~ c(t§ — t e j e )/[(l + z)(l — V cos 9 t )], 9j 
is the half opening angle of the shell, v 1 = T(l — Vcos^)^©, 

S vl oc (i?/i? cro )V~ /3 is the specific spectrum of the radia- 3 The derivatiori of the curvaturc effect has not been presented 

tion in unit of solid angle. On the "equal arriving time sur- in the paper to match the page limit. 



